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Abstract 

Salt, glycerol and dimethyl sulfoxide (DMSO) are used to modify the properties of protein solu- 
tions. We experimentally determined the effect of these additives on the phase behavior of lysozyme 
solutions. Upon the addition of glycerol and DMSO, the fluid-solid transition and the gas-liquid 
coexistence curve (binodal) shift to lower temperatures and the gap between them increases. The 
experimentally observed trends are consistent with our theoretical predictions based on the thermo- 
dynamic perturbation theory (TPT) and the Derjaguin-Landau-Verwey-Overbeek (DLVO) model 
for the lysozyme- lysozyme pair interactions. The values of the parameters describing the interac- 
tions, namely the refractive indices, dielectric constants, Hamaker constant and cut-off length, are 
extracted from literature or are experimentally determined by independent experiments, including 
static light scattering, to determine the second virial coefficient. We observe that both, glycerol 
and DMSO, render the potential more repulsive, while sodium chloride reduces the repulsion. 
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I. INTRODUCTION 



Proteins are complex macromolecules which play a crucial role, e.g. as enzymes or struc- 
tural units, in various processes in living organisms, but also in biotechnology. Many of 
these natural or industrial processes require further ingredients and additives such as salts 
or solvents. In biotechnology, these substances are often added to modify the protein inter- 
actions and phase behavior, for example to favor crystallization [l| . There is a great variety 
of additives; they include ions, such as sodium chloride (NaCl), to regulate the electrostatic 
repulsion or potential, and liquids, such as glycerol or DMSO, to, e.g., stabilize proteins 
against denaturation p], protect proteins against freezing js, 4| or inhibit protein aggrega- 
tion [5]. A thorough understanding of the resulting changes in the interactions and phase 
behavior is beneficial for the design and control of these processes. For example, experience 
has shown that high quality crystals are often formed in slightly supersaturated solutions 
while protein crystals with many defects or amorphous protein precipitates (aggregates) are 
frequently obtained if the sample is quenched deep into the crystal phase, in particular close 
to the metastable gas-liquid coexistence curve where rapid nucleation occurs Thus, 
understanding the phase behavior can guide the design of processes involving proteins. 

Despite the huge complexity of proteins, even of globular proteins such as lysozyme 
considered here, we will use a very simple spherically symmetric model. Although this will 
not do full justice to the details of the protein structure and the individual properties of a 
specific protein, it will allow us to apply concepts developed in soft matter physics to protein 
solutions. It will also shed light on general properties of proteins which can successfully be 
described by a coarse-grained model. On the other hand, the model will fail to describe 
certain aspects and thus indicate what defines the specific properties of a certain protein. 
To describe the protein-protein interactions, we use the Derjaguin-Landau-Verwey-Overbeek 
(DLVO) potential, a well-established model potential in colloid physics 8|, |9(|. Based on this 
interaction model and thermodynamic perturbation theory (TPT) 10|, we will calculate the 
phase behavior of lysozyme in aqueous solution with additives, namely added salt (NaCl), 
glycerol and dimethyl sulfoxide (DMSO), and compare our theoretical predictions to our 
experimental observations. For the additive concentrations studied here, lysozyme does not 



unfold or undergo conformational changes 
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4l3l|. This allows us to test the validity and 



limits of the simple model, to investigate the effects of these additives on the interactions 
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and the phase behavior, and to understand at least quahtatively the trends observed in our 
experiments. 

In several previous studies, the properties of lysozyme solutions, namely the protein- 
protein interactions and the phase behavior, were described using the DLVO potential. Poon 
et al. studied the fluid-solid transition, i.e. the crystallization or solubility boundary, as 
a function of salt and protein concentration as well as the pH of the solution and thus the 
protein charge. They found a universal crystallization boundary if the salt concentration is 
normalized by the square of the protein charge and supported this finding by the behavior of 
the second virial coefficient calculated using the DLVO potential. A more detailed anal ysis 
was later provided by Warren [l5|. Similarly, computer simulations by Pellicane et al. 16 1 
based on the DLVO model showed that the shape of the gas-liquid coexistence curve can 
successfully be predicted, with the critical temperature strongly depending on the Hamaker 
constant. By contrast, Broide et al. 17|] report that the DLVO model cannot fully explain 
their experimentally-observed phase behavior; the predicted height of the Coulomb barrier 
is too low to inhibit protein aggregation. This suggests that hydration forces are important 
to prevent lysozyme from aggregating. 

Despite the widespread use of additives such as glycerol and DMSO, for example in crys- 
tallization essays, only a few studies were concerned with their effect on the phase behavior 
of proteins. Farnum and Zukoski [iS^ studied the effect of glycerol on the protein interac- 
tions in aqueous solutions of bovine pancreatic trypsin inhibitor (BPTI). Upon increasing 
the glycerol concentration, the repulsive interaction increases. A description of these data 
based on the DLVO model requires a Hamaker constant which is much smaller than expected 
from the optical properties of water-glycerol mixtures and of proteins. They thus attributed 
their observations to a small increase of the effective protein size which is caused by an 
enhanced hydration shell (and modeled by an increased cut-off length 6, see Sec. IIII A2p . 
This is supported by observations of Priev et al. 19| . Their experiments indicate an increase 
of the hydration shell upon addition of glycerol and, at the same time, a decrease of the 
core volume of BPTI. Similarly, using small- angle neutron scattering Sinibaldi et al. {20 1 
observed, upon increasing the glycerol content, a significant increase of the hydration shell 
of lysozyme from 0.3 nm to 0.6 nm, and a small decrease of the core volume of lysozyme 



of a 



20 



30ut 6%, such that the total effective protein volume increases by at least about 13% 



2l| . In the presence of glycerol, Esposito et al. 22|] furthermore observed an increase of 
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the hydro dynamic radius with increasing temperature. Using computer simulations, Vage- 



nende et al. 



12| found that glycerol preferentially interacts with the hydrophobic patches 



on the lysozyme surface. Glycerol molecules shield these patches with their hydrocarbon 



backbone, allowing the lysozyme hydration shell to grow. Moreover, Sedgwick et al. [23 1 
studied the effect of up to 40 vol. % glycerol on the phase behavior of lysozyme. They found 
that, based on the DLVO model, the phase behavior can be described almost quantitatively 
if the glycerol-induced changes in the dielectric constant and index of refraction are taken 
into account. The effect of glycerol and DMSO on the fluid-solid boundary of lysozyme 



solutions has also been studied by Lu et al. [2J]. Both additives increase the lysozyme 
solubility. Furthermore, consistent with the previously-mentioned study 23|, Lu et al. 24 1 
found that glycerol and DMSO decrease the protein critical supersaturation, thereby pro- 
moting nucleation. In subsequent work, Gosavi et al. studied the formation of lysozyme 
crystals showing that the crystal growth rate is enhanced upon addition of glycerol. Similar 



observations were^ade by Kulkarni and Zukoski [26 1. Furthermore, Arakawa et al. (3[ and 



Kamiyama et al. [4] observed a preferential hydration of lysozyme in water-DMSO mixtures, 
indicating that DMSO does not bind to the lysozyme surface This suggests that DMSO 
mainly affects the dielectric properties of the bulk solution. Nevertheless, DMSO leads to 
denaturation at volume fractions beyond 70 vol. % [13]. It is hence still not clear whether the 
addition of glycerol and DMSO only changes the solvent dielectric constant and refractive 
index or whether it also has a significant effect on the protein, in particular on its hydration. 

Combining theory and experiment, we investigate the phase behavior of aqueous lysozyme 
solutions with additives. The samples contain 0.7 M or 0.9 M NaCl, up to 20 vol. % glycerol 
and up to 15 vol. % DMSO. By optical inspection, we experimentally determine both the 
fluid-solid transition and the metastable gas-liquid coexistence curves (binodals). Of par- 
ticular interest is their relative location in dependence on the amount of added glycerol 
and DMSO, since protein crystallization is expected to be enhanced in the vicinity of the 
gas-liquid critical point [6]. As we will show, the experimentally-determined phase diagram 
of lysozyme is qualitatively consistent with our theoretical predictions. They are based on 
the DLVO model and thermodynamic perturbation theory (TPT) 10|, and are preformed 
without any free parameters. The values of the parameters describing the DLVO model 
are either taken from the existing literature or determined by independent measurements, 
namely refractometry to measure the refractive indices and static light scattering (SLS) 
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to determine the second virial coefficient. These independent measurements also provide 
detailed information on the effect of the additives on the protein interactions. 

The organization of this paper is as follows. In Section HIl we describe the experimental 
procedures. Section IIII Al introduces the DLVO model for the protein-protein interactions 
and presents the values of its parameters. Based on this model and without additional 
free parameters, in Section IIIIBI the phase behavior is calculated using thermodynamic 
perturbation theory. Then the predicted phase behavior is compared to our experimental 
observations. Finally, our findings are summarized in Section IIVI 

II. MATERIALS AND EXPERIMENTAL METHODS 
A. Materials and Sample Preparation 

We used three-times crystallized, dialyzed and lyophilized hen egg-white lysozyme powder 
(Sigma Aldrich L6876). An aqueous 50 mM sodium acetate buffer was prepared and adjusted 
with hydrochloric acid to pB. = 4.5. The protein powder was dissolved in this buffer to result 
in a concentration of about 40mg/ml. This suspension was several times passed through a 
filter with pore size 0.1 fim (Acrodisc syringe filter, low protein binding. Pall 4611) to remove 
impurities and undissolved protein. This solution subsequently was concentrated by a factor 
of six using an Amicon stirred ultra-filtration cell (Amicon, Millipore, 5121) with an Omega 
10k membrane disc. Pall OM010025. The concentrated solution was used as protein stock 
solution. Its protein concentration was determined in a quartz cell with a path length of 
1 cm by UV absorption spectroscopy at a wavelength of 280 nm using a specific absorption 
coefficient of E = 2.64 ml/ (mg cm). 

Samples were prepared by mixing appropriate amounts of protein stock solution, buffer, 
salt solution (3M NaCl with buffer), glycerol solution (70 vol. % glycerol with buffer), and 
DMSO. Mixing was performed at a temperature above the metastable gas-liquid coexistence 
to prevent phase separation. The glycerol samples contained 0.9 M and 0.7 M NaCl, the 
DMSO samples 0.7 M NaCl, and all samples contained buffer and various concentrations of 
lysozyme. 

In order to obtain protein concentrations above that of the stock solution, we exploited 
the gas-liquid phase separation. A solution was quenched into the metastable gas-liquid 
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coexistence region, centrifuged for 10 min at a temperature below the cloud point temper- 
ature, and the protein-rich phase was collected. The resulting protein concentration of the 
protein-rich phase was calculated from the volume ratio of the two coexisting phases, the 
initial protein concentration and the protein concentration of the protein-poor phase, which 
was determined via its refractive index. 

B. Determination of the Phase Boundaries 

The gas-liquid coexistence curves (binodals) were determined by cloud point measure- 
ments. The sample was filled into a NMR tube with 5 mm diameter, sealed, and placed into 
a water bath at a temperature above the demixing point. The temperature subsequently was 
lowered stepwise and the cloud point temperature identified by the sample becoming turbid. 
This measurement was repeated for samples with different lysozyme volume fractions to 
determine the location of the binodals. 

To determine the crystallization boundary, the sample was filled into an x-ray capillary 
tube and sealed with UV curing glue (Norland). The capillary was then mounted onto 
a home-built temperature stage that allows for observation by optical microscopy (Nikon 
Eclipse 80i). For samples with a low protein concentration, it was necessary to lower the 
temperature to +4°C to induce crystallization. Once crystals were formed, the temperature 
was raised stepwise. The temperature at which the crystals begin to melt was identified with 
the transition temperature. Again, the experiments were repeated with different lysozyme 
volume fractions to obtain the fluid-solid coexistence curves. 

C. Determination of the Refractive Indices 

We used a temperature-controlled Abbe refractometer (Model 60L/R, Bellingham & Stan- 
ley) operated at a wavelength of 589.6 nm to determine the refractive indices of the samples, 
that is of the protein solutions, n, and of the solvent, Ug. The refractive index increments, 
dn/dcp, were obtained from linear flts to the dependence of n on the mass protein concen- 
tration, Cp. 



6 



D. Determination of the Second Virial Coefficient 



Static light scattering (SLS) was performed with a 3D hght scattering instrument (LS- 
Instruments). Due to the low scattering intensity of the samples, the instrument was oper- 
ated with a single beam with wavelength A = 633 nm. The samples were filled into cylindrical 
glass cuvettes with a diameter of 10 mm, centrifuged for at least 10 min at typically 7500g 
prior to the measurements, and placed into the temperature-controlled vat of the instrument 
which was filled with decalin. The samples used for the SLS measurements were dilute with 
lysozyme concentrations between 2.8 mg/ml and 14.4 mg/ml. 

We did not observe a dependence of the scattered intensity, /p, on the scattering angle, 6, 



consistent with the small protein diameter o" = 3.4 nm |27l. 1281 compared to the inverse of the 
scattering vector q = {Ann/ \) sin (9/ 2) so that qa <t^l |29|. Therefore, the time-averaged 
intensity scattered by the protein solution, (/p(cp)), was detected at a single scattering angle, 
6 = 90°, corresponding to a magnitude of the scattering vector q ^ 0.018 nm~^. The excess 
scattering due to the protein is expressed as the Rayleigh ratio. 



i?(cp) ^ , (1) 



which relates the difference in the average scattering intensities of the protein solution and 
the solvent, (/p(cp)) — (Jg), to the scattered intensity of the toluene reference, {It), with 
known Rayleight ratio, Rt- The ratio of the refractive indices of the protein solution, n, 
and toluene, rit, accounts for the different sizes of the scattering volumes in the two cases. 
The temperature dependence of Rt was experimentally determined from the temperature 
dependence of the intensity scattered by a toluene sample. The absolute value Rt = 1.40 x 
10~^cm^^ at T = 35°C was taken from literature [3^ and is consistent with other reports 
(see |31|] and references therein). 

For dilute solutions, R{cp) is related to the second virial coefficient, B2, by 

Kcp 1 f_ , 2NaB2 



l + ^T^Cph (2) 



i?(cp) M(o) V MW 

where M*^*^) and M^^^ both represent the molar mass M, but by two different values (see 
below), A^A is Avogadro's constant and K an optical constant given by 



Note that the refractive indices and the refractive index increments in the above equations 
depend on temperature. 

SLS measurements of a series of protein concentrations hence allow to determine B2 
from the slope of the Cp dependence, and the molar mass M from the intercept at Cp — 0. 
The measurements were performed at different temperatures T to obtain the temperature 
dependence of B2(T), while M is expected to be temperature independent. The values of 
the molar mass, experimentally determined at different temperatures from the intercept at 
Cp — )■ 0, lie in the range 14 800 g/mol < M'^'^^ < 19 500 g/mol, while the expected value is 
14 320g/mol Q- The error in M^°) is mainly due to the low scattered intensity, uncertainties 
in dn/dcp and i?f Dynamic light scattering experiments on the same dilute samples did not 
indicate the presence of aggregates, which are also neither expected for aqueous solutions 



E 



33| nor for glycerol-water 



12| or DMSO-water mixtures 13|. To avoid that experimental 
errors in the molar mass affect B2{T), we used the literature value 14 320 g/mol in the term 
for the slope, 2NaB2/M^^\ i.e. M^^^ = 14 320 g/mol, when calculating B2{T) based on 
Eq. ([2]) [34]. This renders B2{T) independent of the absolute scattered intensity, i.e. M*^°\ 
which has a significant uncertainty, and thus makes B2{T) more reliable. 

III. RESULTS AND DISCUSSION 

We first present the simple model used to quantify the protein-protein interactions, 
namely the DLVO potential model (Sec. IIII A 1|) . as well as the values of its parameters. 
These are either taken from the literature, or determined by independent measurements 
(Sec. IIII A 2|) . Based on this interaction potential and these parameter values, we calculate 



the phase diagram without any free parameter (Sec. IIII B ip . Our experimentally determined 
phase diagrams (Sec. IIII B 2h are then compared to the theoretical predictions (Sec. IIII B 3h . 
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A. Protein-Protein Interaction Potential 



1. DLVO Model 

The DLVO pair potential, u{r), includes three parts; hard-sphere UQ{r), repulsive elec- 
trostatic Uei{r) and attractive van der Waals interactions u^dwi^), 

u{r) = Uo{r) + Mei(r) + Mvdw('^) , (4) 

where r is the center-to-center distance of two spherical particles. 
The hard-sphere contribution is 

Mr) = { (5) 
I 0, r>a. 

The electrostatic repulsion is described by a screened Coulomb potential 

/3Mei(r) = —2 > (t) , (6 

{l + Ka/2) r 

where /3 = l/{kBT) with fee the Boltzmann constant, cr the diameter, Z the charge and Ib = 
I (47r£:o£^s^B7') the solvent-specific Bjerrum length depending on the elementary charge, e, 
dielectric constant in vacuo, eo, and dielectric constant of the solvent, Eg. The square of the 
electrostatic (Debye) screening parameter, k, is given by 

where = Cp/ pp is the particle volume fraction with pp the mass density, Cg the molar salt 
concentration and Cb the molar concentration of dissociated buffer. The factor 1/(1 — 0) 
corrects for the volume occupied by proteins and thus takes into account the free volume 



accessible to the micro-ions, which cannot penetrate the protein [35|, l36|. We ignore here 
the small free volume corrections arising from the finite size of the microions. 
The van der Waals attraction is of the form ^ 

.2" 



Mvdw(r) = -— ^ + — + 2 In 



{r>a + S), (8) 



where A is the (effective) Hamaker constant and 6 a cut-off length. The cut-off length S 
represents the smallest possible surface-to-surface separation of two proteins, introduced to 
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avoid the unphysical l/{r — a) divergence at r — )■ cr. Its value can be estimated by the 



thickness of the Stern layer, i.e., by the mean size o: 
The Hamaker constant can be approximated by 



;he salt ions. 



A = A^=o + A,,>o = -:kBT 



+ 



[nz 



2^2 



16v^ (nl + nif^ ' 



(9) 



where Ay=Q is the zero-frequency and A^^q the frequency- dependent part, h the Planck 
constant, and the static dielectric constants of the solvent and protein, and 
their refractive indices, and z/g = 3 x 10^^ the characteristic UV electronic adsorption 
frequency js]. 



2. Values of the Model Parameters 



To quantitatively determine the potential M(r), we require the values of the refractive 
indices and static dielectric constants, of the cut-off length 5, and of the parameters describ- 



ing the protein 



cr = 3.4 nm 



27 



ysozyme is an approximately spherical protein with an effective diameter 



28(1, molar mass M = 14 320g/mol, mass density pp = 1.351 g/cm'^ 32| 
and Z = 11.4 net positive elementary charges at the present pR = 4.5 [s^]. The acidity 
constant of sodium acetate is pK^ = 4.76 resulting in a dissociated buffer concentration of 
Cb = 27 mM at pH = 4.5. 

a. Refractive indices and static dielectric constants. The protein index of refraction is 
obtained by the linear extrapolation. Hp = Q M / {Np^Tia^) dri/dcp + ns, where 6/{Natt(J^) is 
the molar concentration of a 'bulk' protein. We find no significant dependence of dn/dcp 
on composition or temperature in the relevant parameter range, consistent with previous 



findings 
1.72. 
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Averaging dn/dCp over all tested temperatures and compositions yields rip 



It is usually assumed that dry protein has a static dielectric constant in the range of 



2 < £p < 4 



23 



391]. However, measurements with protein powders often indicate larger 



values, for example 1 < < 10 for wet lysozyme powder [40[, 20 < ^p < 37 for hydrated 



lysozyme crystals 



4l| and, based on MD simulations, ^p = 25.7 for lysozyme in aqueous 



solution [42] . The larger values are attributed to water or moisture leading to hydration and 



some water penetration 



39|. For consistency with previous work 23|, |39|], we use a static 
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FIG. 1: Temperature dependence of the index of refraction of the solvent, ns(T), for different added 
NaCl concentrations, Cg, as well as glycerol and DMSO contents (as indicated). The lines are linear 
fits to the experimental data points. 

dielectric constant of lysozyme given by = 2. We might underestimate in this way and 
thus overestimate the static contribution to the Hamaker constant Ai,=q (Eq. [9]) and hence 
the protein attraction. However, the effects are not significant, because, as will be shown 
below, the main contribution to the Hamaker constant is due to the non-zero frequency term 

The refractive index of the solvent, Ug, was measured and found to decrease with increas- 
ing temperature (Fig. [1]). The refractive indices of glycerol, rig, and DMSO, noMSO; are 
larger than the one of water, n^. For example, rig = 1.47, wdmso = 1-48 and = 1.33 at 
T = 20°C. The addition of glycerol or DMSO thus increases rig. 



The dielectric constant of the solvent, e,, was taken from literature 
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49j (Fig. ED. It 



decreases with increasing temperature, which is due to the increasing thermal motion of the 



11 




FIG. 2: Temperature dependence of the die! 
(a) water-glycerol |43l445l| . (b) water-DMSO 



ectric constant, £s{T), of different solvent mixtures. 



47 



491 ]. The hues are 



4g], (c) aqueous salt solutions 
linear fits. The fits that belong to solutions containing 10 and 20 vol. % glycerol are based on five 
data points in an extended temperature range 20°C <T< 80°C. 

solvent molecules and reduces the alignment of the dipolar molecules by an external electric 
field. The dielectric constants of glycerol, Eg, and DMSO, eoMso, are much smaller than the 
one of water, for example eg = 41, ^dmso = 48 and = 80 at T = 20°C. Therefore, the 
solvent dielectric constant, ^s, decreases with increasing glycerol and DMSO content, with 
the decrease being more pronounced for glycerol (Figs. and b). There is also a significant 
decrease of with increasing salt concentration, Cg] adding 0.7 M NaCl lowers Eg more than 
adding 20 vol. % glycerol (Fig. [2t) 47H49| . This decrease is related to the decrease in the 
electric polarizability of the solution caused by the alignment and attraction of the water 
molecules by the salt ions. 

Based on these static dielectric constants, we calculate the repulsive electrostatic inter- 
action part, Uei{r), (Eqs. and ([7])) for the different compositions, i.e. different added 
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FIG. 3: Temperature dependence of the contact value of the reduced electrostatic interaction 
potential, /3nci(cr) = Z'^lB/{cr[l + ko"/2]^), for different added NaCl concentrations, Cg, as well as 
glycerol and DMSO volume fractions (in vol.%, as indicated). Note the break in the ordinate axis. 

salt, glycerol and DMSO content. Upon increasing the salt (NaCl) concentration Cg, the 
screening parameter, k, increases and thus the contact value of the electrostatic potential, 
Mei(cr), decreases (Fig. [3]). In addition, Uci{<7) depends on the solution composition since the 
Bjerrum length /b depends on Es- Adding glycerol or DMSO decreases Eg (Figs. [2ti and b) 
and hence increases Uei(cr), with the effect of glycerol being significantly stronger. 

We have determined the zero- and non-zero frequency parts of the Hamaker constant 
(Eq. iQ) based on the above values of the refractive indices and static dielectric constants. 
The non-zero frequency part, A^^q, is about ten times larger than the zero-frequency part, 
Ai,=o, and hence dominates A (Fig. H]). Both parts decrease with increasing temperature as 
well as increasing salt, glycerol and DMSO content. Glycerol and DMSO cause about the 
same decrease of Ay>o, reflecting a similar effect on Us (Eq. ([9]) and Fig. [1]). In contrast, we 
find a less significant decrease of A,y=Q upon addition of DMSO than of glycerol, due to the 
weak effect of DMSO on Eg (Fig. [2}d). Added salt (NaCl) causes a pronounced decrease of 
Ai,=o since salt significantly lowers Es (Fig. [Jt). For our sample compositions, i.e. solvent 
mixtures with added salt, the temperature dependence of Es is not available. However, since 
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FIG. 4: Temperature dependence of the (a) zero-frequency, Ajy=o, and (b) non-zero frequency, 
^u>o, parts of the Hamaker constant in Eq. ([9]) for different compositions, (a) water-glycerol 
mixtures without NaCl (purple dashed hnes), water-DMSO mixtures without NaCl (green sohd 
hues), aqueous salt solutions with 0.7 M NaCl (black dotted line) and 0.9 M NaCl (blue solid line), 
(b) water-glycerol mixtures containing 0.7 M NaCl (black dotted lines), water-glycerol mixtures 
containing 0.9 M NaCl (blue solid lines), water-DMSO mixtures containing 0.7 M NaCl (green 
solid lines). The numbers indicate the volume fractions (in vol.%) of glycerol and DMSO. 

Ai,=o is much smaller than A,y^Q, we neglect the effect of added salt and consider only the 
effect of glycerol and DMSO. Nevertheless, in our calculations we take into account the effect 
of added salt on Ug and thus A^^q (Fig. HJd), which represents the dominating contribution 
to the Hamaker constant. Glycerol and DMSO not only similarly affect A^^q and thus the 
Hamaker constant A, but also the total DLVO pair potential, u{r), (see Fig. [5]). 

b. Cut-off' length 6. There is one potential parameter left to be determined, namely the 
cut-off length, 6. We choose it such that it reproduces the second virial coefficient, B2{T), 
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FIG. 5: DLVO pair potential, u{r), for NaCl concentration Cg = 0.7 M, protein volume fraction 
<j) = 0.002, and temperature T = 20° C, for an aqueous solution (black line), a water- glycerol 
mixture (lOvol.%, red line), and a water-DMSO mixture (5vol.%, green line). The vertical lines 
indicate the cut-off distances 6 of the van der Waals attraction. 

which we have experimentally determined by static light scattering. 

From a theoretical point of view, the second virial coefficient, B2{T), is introduced by 
expanding the reduced free energy density, /(T, 0), in a power series of the particle volume 
fraction, 0, according to 

/(T, <f>) = f,^{<jy) + 4h{T)<jy' + 0{<jy') , (10) 

which involves the ideal gas part of the free energy density, Jq'^, and its leading-order cor- 
rection, the normalized or reduced second virial coefficient 

Mr) = ^. (11) 

The latter is equal to the ratio of the second virial coefficient, 

B2 = -^jdr (exp [-(3u{r)] - 1) , (12) 

to the second virial coefficient of a hard-sphere system, B^^ = 27r((T + S)^/3. For the DLVO 
pair potential, u{r), (in Eq. (jl])) the normalized second virial coefficient, b2, can be written 

15 




FIG. 6: Temperature dependence of the normalized second virial coefficient, 62 (T), for different 
added NaCl concentrations, Cg, as well as glycerol and DMSO contents (as indicated), experi- 
mentally determined by static light scattering. The solid lines represent fits assuming a linear 
temperature dependence of the cut-off length 5. 



as 

"OO 



3 /"°° 

^2 = 1- / ''^ (e^P [-/5(«ei(r) + Mvdw(r))] - 1) . (13) 

The second virial coefficient characterizes the thermodynamic strength of the pair interac- 
tions. Positive values of 62 indicate that repulsive interactions dominate, while negative 
values correspond to dominating attraction. 

Using SLS, we determined the temperature dependence of 62 (r) in the range 10°C < 
T < 40°C for all relevant solution compositions, i.e. salt (NaCl) concentrations, Cg, as well 
as glycerol and DMSO volume fractions (Fig. [6]). Outside this temperature range, reliable 
SLS experiments cannot be performed. With increasing temperature, as expected h2{T) and 
thus the repulsion increases, which reflects the increasing electrostatic (Fig. [3]) and decreasing 
van der Waals interactions (Fig. Hj). This leads also to an increase of h2{T) with increasing 
glycerol and DMSO concentrations, in agreement with previous observations 23|, |26|, |50|, 



with DMSO having a signiflcantly more pronounced effect. Furthermore, h2{T) increases on 
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lowering Cs, due to the increased electrostatic repulsion (Figs. [6^ and b). The sample with 
Cs = 0.7 M and 20 vol. % glycerol even shows a positive b2(T) in the investigated temperature 
range 10°C < T < 40°C, indicating that the lysozyme pair potential is predominantly 
repulsive at high glycerol concentrations. 

The experimental b2{T) is fitted by the DLVO potential prediction (Eq. (IT^ ). with the 
cut-off length 6 being the only adjustable parameter. A constant 6 does not reproduce the 
experimental data satisfactorily. We thus assume a linear temperature dependence of the 
cut-off length, 6(T) = 60 + 61T, and fit the coefficients 60 and 61 to obtain optimal agreement 
between experimental and theoretical b2{T) in the temperature range 10°C < T < 40°C 
(Fig. E]). The fitting is done by applying the Levenberg-Marquardt method 5l(|, and the 
improper integrals appearing in the expression for the second virial coefficient are evaluated 
using the Chebyshev quadrature. For all considered compositions and the temperature range 
of the light scattering experiments, S{T) increases with increasing T (Fig. [7]). Its values are 
in the range 0.15 nm < S{T) < 0.3 nm, except for the solution with 0.7 M NaCl and 20 vol. % 
glycerol for which unphysically large S{T) are obtained. This range of values agrees wit 
previously published values, which are in the range of 0.1 nm < S{T) < 0.3 nm 
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28 



52 



54| . The large values for 6{T) obtained for 0.7 M NaCl and 20 vol. % glycerol are due to 



the positive- valued experimental ^2(^)7 while the DLVO model predicts a predominantly 
attractive potential and thus a negative h2{T). This forces 5(T) to rapidly grow with T to 
reduce the attractive van der Waals contribution. The effect of the van der Waals interaction 
on h2{T) is controlled by the Hamaker constant A and the cut-off length 5. A change in one 
of the two hence can be compensated (or enhanced) by a corresponding change in the other 
parameter. 

The cut-off length, 5, accounts for molecular effects on a coarse-grained level, such as the 
solvation shell, which might contain both water and additive molecules (see Sec. [T]), or the 
condensation of counterions on the protein surface, i.e. in the Stern layer. Its increase with 



12 



20 



m. 



increasing temperature and glycerol content is consistent with previous findings 
As mentioned above, glycerol affects h2{T) less than DMSO, i.e. the addition of 10 vol. % 
glycerol has about the same effect as 5vol.% DMSO (Figs. |6]3 and c), although the optical 
properties of the two solvents suggest a similar effect on A^^q and thus A (Fig. H]) and, 
due to the small effect of glycerol and DMSO on the electrostatic interactions (Fig. [3]), also 
on the DLVO potential (Fig. [S]). This results in a larger cut-off length 6{T) for DMSO 
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FIG. 7: Temperature dependence of the cut-off length S(T) as obtained from fits of DLVO-based 
predictions to the measured 62 (T). (a) Different added NaCl concentrations, Cg, and glycerol 
contents and (b) same NaCl concentration Cg = 0.7 M, but different glycerol and DMSO contents 
(as indicated). 

for a given concentration of glycerol or DMSO (Fig. [7]), and suggests that glycerol and 
DMSO interact differently with lysozyme. Compared to glycerol, little is known about the 
DMSO-lysozyme interactions, apart from a preferential hydration of lysozyme in water- 
DMSO mixtures indicating a weak or no binding of DMSO to lysozyme [3|, |4j. This suggests 
that DMSO mainly influences the dielectric properties of the solution, and is consistent 
with the observed weaker temperature dependence of the cut-off length S{T) in the case of 
DMSO. 

B. Phase Behavior 

We have now determined all values of the parameters determining the DLVO model. We 
can thus predict the phase behavior without any free parameters and compare it to our 
experimental observations. 
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1. Phase Behavior and Free Energy by Thermodynamic Perturbation Theory 



Our prediction of the phase behavior of lysozyme is based on expressions for the Helmholtz 
free energy which are obtained by the second-order thermodynamic perturbation theory 
(TPT) of Barker and Henderson lOj. Knowledge of the free energy of the fluid and solid 
phases allows us to deduce all equilibrium thermodynamic properties, such as the fluid-solid 
and gas-liquid coexistence curves including the critical point. The Barker-Henderson theory 
is known to yield quite accurate results IsS, IsG]. It is therefore widely used to calculate the 



phase diagram of protein solutions 



57|-|61| 



Within TPT, the free energy, F, is expanded in powers of the perturbation potential 
Up = u — Uq, where Uq is the pair potential of a reference system and u the total protein- 
protein pair potential (see Eq. (jl])). The reference system is an effective hard-sphere system 
with diameter a + 6, which takes into account the cut-off length S in Eq. ([8]). Instead of the 
free energy, F{N, V, T), of a system with proteins in a volume V, we use the dimensionless 
free energy density / = (SFvq/V, where Vq = na^/G is the particle volume. Expanding the 
free energy density, /, up to second order in Up gives 

12 



/(0,T) = /o(0) + 



(a + 5y 
6 



0' / drr go{r) /3up{r) 



dr go{r)[f3up{r)y 



(14) 



where 0' = 0(1+5/0")^ is the effective sphere volume fraction, and XT/(/3'i'o) = l/0(<90/5Ho)rp 
the isothermal compressibility of the reference system with Hq its dimensionless osmotic 
pressure. Furthermore, go{r) is the radial distribution function of hard spheres in the fluid 
phase, and the orientationally-averaged pair distribution function of hard spheres in the 
solid phase. For the radial distribution function in the fluid phase, we use the Verlet- 



Weis corrected 1621 Percus-Yevick solution 
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6J], and for the orientationally-averaged 



pair distribution function the expression by Kincaid 65|] for a face-centered cubic (fee) 
crystal phase. The free energy density, /o(0), of the effective hard-sphere reference system 
is separately defined in the fluid and solid branches. In the fluid phase, /o(0) consists of two 
parts: First, the ideal gas part. 



«/^'(0) = 0' [ln(0'AVt^o) - 1] 



(15) 



where A = /i/y 27rmp/cBT is the thermal wavelength, nip the protein mass and a = 0'/0. 
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Second, the interaction part which is approximated by the accurate Carnahan-Starhng equa- 
tion of state [66,] 

In the sohd phase, fo{4>) is described by Wood's equation of state \^ on assuming a fee 
crystalhne lattice 

atrW ^ 2.1306^' H- 3^'.„ (^-|^) + (^) . (17) 

where 0cp = 7r-\/2/6 is the fee volume fraction for closed packing. To evaluate the free energy 
density, f{(f), T), in the solid phase, the integrals in Eq. f|T^ are divided into intervals centered 
around the crystal lattice sites, which is necessary to sufficiently resolve the discreteness of 
go{r) that develops with increasing 0. 

The coexistence of gas (g) and liquid (1) phases requires an identical temperature T, 
osmotic pressure 11 and chemical potential fi in both phases. While the temperatures of two 
phases in contact are always identical under equilibrium conditions, the latter two conditions 
determine the two volume fractions, namely that of the coexisting gas, 0g, and liquid phase, 



ng(T, 0g) = ni(T, 00 with n(T, 0) = 0^ (^^^^^^^|^) ^ (is) 

/ig(T, 0g) = /ii(T, 0i) with /3/i(T, 0) = (^^^^^) ^ ■ (19) 
Similarly, the coexistence of fluid (f ) and solid (s) phases is determined by the two conditions 

nf(T,0f) = n,(T,0s) (20) 
/if(r,0f) = ^s(T,0s), (21) 

where 0/ is the volume fraction of the fluid, and 0^ the one of the coexisting solid phase. 
To compute the phase coexistence curves, the Newton-Raphson method with line search is 



applied [51 1. 



The critical point terminates the gas-liquid coexistence curve. It is determined by the 
vanishing second and third derivatives of the free energy density, i.e., by 

and ^^0. (22) 
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The higher-order derivatives of the free energy are calculated using Ridder's implementation 



of Neville's algorithm 



5l|. 



Before comparing our theoretical and experimental findings, we test the accuracy of our 
TPT implementation. No simulated phase diagrams are available for the DLVO model 



with our parameter values. However, Monte Carlo simulations 55|, |68|] have been used to 



determine the phase diagram of hard spheres interacting via an attractive Yukawa potential 

u{r) ={ 23 

where ^ > is the (contact) strength, and the range of the attraction. 

Using TPT, we calculate the phase diagram of this system. The calculated fluid-solid co- 
existence curves agree with the simulated ones for all investigated values of (see Fig. |8]). 
Also the widths of the gas-liquid coexistence curves are well reproduced. The calculated 
gas-liquid binodals, however, are located slightly above the simulation data. We suppose 
that the binodals are shifted because the mean-field type TPT does not account for critical 
density fluctuations. Critical fluctuations allow a system to probe density inhomogeneities 
and thus tend to stabilize the homogeneous supercritical fluid phase against phase separa- 
tion, therefore shifting the critical point to lower temperatures. A renormalization-group 



correction scheme has been derived 



69| which, according to previous calculations [70|, is 



expected to flatten the coexistence curves around the critical point. 



2. Experimentally Determined Phase Behavior 

We experimentally determined the fluid-solid and gas-liquid coexistence curves for dif- 
ferent contents of added salt (NaCl), glycerol and DMSO (Fig. [9]). For all compositions 
explored, the gas-liquid binodals lie within the fluid-solid coexistence regions and are thus 



metastable. This is typical for short-ranged attractive pair potentials IT], I7ll-l73| . The bin 



odal widens rapidly upon cooling at volume fractions (f) below the critical volume fraction 
(pci while the data for (p > (pc suggest almost flat binodals at volume fractions larger than the 
critical point value. Similar shapes of lysozyme binodals were observed previously {3, S|- 
The effect of salt on the phase behavior was studied for 0.7 and 0.9 M NaCl (Figs. |9K 
and b, respectively). With increasing salt concentration, Cg, the fluid-solid and gas-liquid 
coexistence curves move to higher temperatures, and the width of the gas-liquid binodal 
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FIG. 8: Phase diagrams with fluid-soUd and gas-hquid coexistence curves of attractive Yukawa 
potential systems with different ranges of attraction, z"^, as indicated. Coexistence curves with 



ack filled circles, and the 
6811 ■ The coexistence 
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critical points calculated using TPT are represented as solid lines with b 
data from Monte Carlo simulations are represented as red open circles 
curves are shown in the normalized density (/oo"^) - normalized potential strength (/3^) - plane. 
Note that the (/3^) - axes run downwards. Furthermore, no simulation data for the fluid-solid and 



gas-liquid coexistence curves for z = 5.0 and z = 9.0, respectively, were provided in 



55 



increases. These findings are consistent with the electrostatic screening by added salt. Upon 
increasing Cg, the screening length, k^^, decreases. Thus the electrostatic repulsion decreases 
and attraction becomes more important. This favors crystals compared to the fluid phase 
and also the (metastable) gas-liquid phase separation compared to a homogeneous fluid. The 
fluid-solid as well as the gas-liquid coexistence curves hence shift to higher temperatures, 
and the gas-liquid coexistence curve becomes broader, as experimentally observed. 

Upon addition of glycerol and DMSO, both coexistence curves shift to lower temperatures 
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FIG. 9: Experimental phase diagram of lysozyme solutions with different added NaCl concentra- 
tions, Cs, as well as different glycerol and DMSO contents (as indicated), in dependence on the 
volume fraction, (p, and temperature, T. The open circles, inverted triangles, pentagons, stars and 
rotated triangles indicate the fluid-solid coexistence curves. The open squares, triangles, diamonds, 
hexagons and octagons indicate the gas-liquid coexistence curves. The filled squares, triangles and 
diamonds indicate the critical points. 



(see Fig. [H]). The shift of the binodals is more pronounced, which widens the gap between 
the two coexistence curves. This trend has been observed before in other systems and is 
characteristic for a decreasing range of attraction 



55|,[3l (Fig. [5]). The critical point is not 
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only shifting to lower temperatures, but also to slightly lower volume fractions. Considering 
the quantitative effect on both curves, we find that the addition of 10 vol. % glycerol and 
5vol.% DMSO induce a similar shift. The higher efficiency of DMSO is consistent with its 
more pronounced effect on b2{T) (Fig. [6]). 

We have investigated the relation between the second virial coefficient, b2{T), and the 
location of the fluid-solid and especially the gas-liquid coexistence curves. For this purpose, 
we parameterized the temperature dependencies of the experimentally-determined second 
virial coefficients (Fig. [6]) by the phenomenological relation b2(T) = ao + «! T + a2 T^, where 
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0.9, 20% Glyc. 
0.7, 20% Glyc. 
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FIG. 10: Experimental phase diagram of lysozyme solutions with different added NaCl concentra- 
tions, Cg, as well as different glycerol and DMSO contents (as indicated), as a function of volume 
fraction, (j), and reduced second virial coefficient, 62- The second virial coefficient has been (a) ex- 
perimentally determined and (b) calculated based on the DLVO model. The filled symbols indicate 
the fluid-solid (fs) coexistence curves, the open symbols the gas-liquid (gl) coexistence curves. 



.eads to a collapse of 
Q,y,Q- This 



ai and 0:2 depend on the composition of the solution. Using this relation, the binodals are 
plotted in dependence of 62 and in Fig. [TOk . This representation 
all binodals onto a single curve, as it has been observed previously 
indicates that in our system the binodals are controlled by 62, i-e. by a single integral 
parameter characterizing the potential, and not by the details of the potential, such as its 
depth and range. By contrast, this scaling of 62 is not observed for the fluid-solid transition. 



Nevertheless, this transition occurs for 62 values consistent with previous reports [77 1. 
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If instead the 62 values calculated based on the DLVO model are used, the fluid-solid and 
gas-liquid coexistence curves do not scale (see Fig. [TOb) and the critical points are different. 
This indicates that the scaling very sensitively depends on the temperature dependence 
of b2{T) and thus on the model. Small errors might be introduced by the simpliflcations 
implicit in the DLVO model. However, the main effect is due to the extrapolation of b2{T) 
beyond the available temperature range 10°C < T < 40°C (see Fig. [6]), which indicates the 
need for a very reliable and robust parameterization of the temperature dependence of the 
cut-off length S{T). To avoid spurious effects, in the following we scale the temperature and 
concentration axes by the critical temperatures and concentrations. This axis scaling also 
reduces artifacts due to the slight shift in the binodals caused by the mean-field nature of 
TPT (Fig.El). 

3. Comparison of Calculated and Experimentally Determined Phase Behavior 

The experimental phase diagrams are now compared to our TPT predictions based on 
the DLVO model. The axes are scaled by the critical temperature, Tc, and volume fraction, 
0c, which matches all theoretical and experimental critical points (see Fig. [TT]) . (Note that 
in this representation the fluid-solid transition shifts to higher T/T^, in contrast to the 
representation using the absolute temperature scale (Fig. [H]). This representation reveals 
an increase of the gap between the fluid-solid and the gas-liquid coexistence curves upon 
addition of glycerol and DMSO for the experimental data and, less pronounced, in the range 
<f)/4>c < 0.2 for the theoretical predictions. Note that in the main part of the phase diagram 
the theoretically predicted gap actually decreases. In view of the successful test of the TPT 
approach (Fig. [8]), we attribute this discrepancy to the inherent approximations of the DLVO 
potential and in particular to the model for the cut-off length 6. 

The width of the gas-liquid coexistence curve is experimentally found to remain practi- 
cally unchanged in this representation when glycerol or DMSO is added, while theory pre- 
dicts a slight increase of the width. The theoretically and experimentally observed widths 
quantitatively agree at the largest glycerol content of 20 vol. % (Fig. [11]). In particular, our 
calculations do not reproduce the flat binodal region around the critical point, which we 
observe experimentally. As discussed above (see Fig. [H]), this is attributed to the mean-fleld 
character of the theory. This is furthermore consistent with previous experiments [3] and 
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FIG. 11: Experimental and predicted phase diagram of lysozyme solutions with different added 
NaCl concentrations, Cg, as well as different glycerol and DMSO contents (as indicated), as a func- 
tion of reduced volume fraction, 0/</>c, and reduced temperature, T/Tc, where {(pcTc) defines the 
critical point. Solid lines represent the TPT calculations based on the DLVO model, open sym- 
bols the experimental data for the stable fluid-solid (fs) and metastable gas-liquid (gl) coexistence 
curves. Filled symbols mark the critical point. 



simulations [16i]. In particular, similar gas-liquid coexistence curves have been theoretically 
observed for colloidal dispersions with competing repulsive and attractive interactions, and 
have been attributed to critical fluctuations which become significant in an enlarged region 



of the phase diagram 78 1. 
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IV. CONCLUSIONS 



In a combined theoretical and experimental effort, we investigated the phase behavior of 
aqueous lysozyme solutions in the presence of additives, namely glycerol and DMSO. We in 
particular investigated the fluid-solid and gas-liquid coexistence curves. Upon the addition 
of glycerol or DMSO, both curves shift to lower temperatures and, in addition, the gap 
between these two curves increases, consistent with a decreasing range of attraction upon 
increasing the content of additives. 

Our theoretical calculations are based on a DLVO interaction potential. This models the 
proteins as hard spheres with an isotropic electrostatic repulsion, which is caused, e.g., by 
uniformly-distributed charges, and an isotropic van der Waals attraction characterized by 
the (effective) Hamaker constant and a cut-off length. The Hamaker constant is calculated 
based on the experimentally determined, or in the literature available, macroscopic optical 
and dielectric properties of the solvent mixtures and the protein, namely their indices of 
refraction and static dielectric constants. The cut-off length 6 was adjusted to reproduce the 
second virial coefficient, which was experimentally determined by static light scattering. We 
find cut-off lengths in the range 0.15 nm < S{T) < 0.3 nm, with the value increasing with 
increasing temperature and glycerol content and, even more pronounced, with increasing 



21, 
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52N54| . The cut-off length subsumes 



DMSO content, consistent with previous studies 
molecular effects, such as a hydration shell or condensed counterions on the protein surface. 
Our experiments thus indicate that glycerol and DMSO interact differently with lysozyme. 
DMSO seems to affect the lysozyme interactions by mainly lowering the solvent dielectric 
constant, while glycerol likely changes the hydrophobicity of the lysozyme surface. This is 



12|. 



in accordance with recent computer simulation results 

Having determined the values of all parameters, we calculated the interaction potential 
for different temperatures and compositions, added salt (NaCl) concentrations as well as 
different glycerol and DMSO contents. We found that glycerol and DMSO increase the 
repulsive interaction, as quantified by the second virial coefficient 62- For a sufficiently high 
content of additives, a predominantly repulsive protein interaction is observed. The increase 
of the repulsion is significantly stronger for DMSO, although the optical and dielectric prop- 
erties of glycerol and DMSO are very similar. We hence attribute this difference to the 
specific molecular interactions of the additives with lysozyme, characterized by the cut-off 



27 



length S{T). 

Based on the DLVO model and its parameters, we predict the phase behavior without 
any free parameters. Our calculations are based on the thermodynamic perturbation theory 
(TPT), which was successfully tested and found to be numerically fast. The theoretical 
predictions were then compared to our experimental observations, for the fluid-solid and 
the gas-fiuid coexistence curves. Particular attention was paid to the effect of additives on 
these curves, and the size of the gap between them. We find that the general trends in the 
phase behavior, and especially the effects of the additives, are also reproduced by the DLVO 
model, although not on a quantitative level. 

The strength of the DLVO model lies in its simplicity, for it allows to incorporate the effect 
of solvent mixtures by taking into account their macroscopic optical and dielectric properties. 
More complex models are available, but might need to be adapted to proteins and solvent 
mixtures or additives. For example, theoretical calculations based on anisotropic attractive 
interactions successfully describe both the width of the binodal and the gap between the 
binodal, and the fluid-solid transition for aqueous protein solutions without additives 60, 



79 



80| . Anisotropic attractive interactions are indeed induced by hydrophobic patches on 



the lysozyme surface, and have been taken into account in a patchy sphere model by Curtis 
et al. |8l|. Anisotropic attractions as well as repulsions, as due e.g. to an anisotropic 



charge distribution or van der Waals interaction, are not taken into account in our DLVO 
model. Our model furthermore neglects non-electrostatic contributions to the lysozyme pair 
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59|. 



potential that originate from dispersion forces between the micro-ions and proteins 
The range and strength of such non-electrostatic and hydrophobic interactions has not been 
determined experimentally, and hence cannot be taken into account unambiguously, e.g., 
through an effective potential of mean force. Furthermore, hydrophobic interactions, as well 
as non-rigid and non-spherical shapes, specific surface properties, ion-specific interactions, 
hydration, and other effects need to be explored in future studies. 
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